Draft version February 2, 2008 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



MOLECULAR GAS AND NUCLEAR ACTIVITY IN ULTRALUMINOUS INFRARED GALAXIES 

WITH DOUBLE NUCLEI 

A. S. Evans 1 ' 2 , J. M. Mazzarella 3 , J. A. Surace 4 , & D. B. Sanders 5,6 

Draft version February 2, 2008 



(N 
O 
O 
(N 

bJQ' 
< 

as 



in 
oo 
o 

O 

6 

c3 



ABSTRACT 

High-resolution C0(1 — > 0) observations of five ultraluminous infrared galaxies (ULIGs: Lm[8 — 
1000/im] ;> 10 12 L Q ) with double nuclei arc analyzed. These sources constitute a complete subset of 
local ULIGs expected to be in an intermediate stage of merging and selected with projected nuclear 
separations of 2'.'0 — 5'.'4 (3-5 kpc) so they could be resolved with the Owens Valley Radio Observatory 
Millimeter Array. The observed pairs include two mergers with cool far-infrared colors (25/im to 60/im 
flux density ratio /25yum//60/im < 0-2) from the Infrared Astronomical Satellite (IRAS) Bright Galaxy 
Sample (IRAS 12112+0305 and IRAS 14348-1447) and three mergers with warm infrared dust temper- 
atures (/25pm //eopm > 0.2) selected from the IRAS Warm Galaxy Sample (IRAS 08572+3915, IRAS 
13451+1232 = PKS 1345+12, and IRAS 13536+1836 = Mrk 463). These ULIGs are further distin- 
guished by the presence of pairs of active nuclei; among the ten nuclei, nine have Seyfert or LINER 
classifications and one is unclassified. 

Molecular gas is detected only on the redder, more radio-luminous nucleus of the warm objects, whereas 
both nuclei of the cool ULIGs are detected in CO. The inferred molecular gas masses for the detected 
nuclei are 0.1 — 1.2 x 10 10 M Q , and the undetected nuclei have molecular gas masses at least 1.2-2.8 times 
less than that of their CO-luminous companions. Upper limits on the extent of the CO emitting regions 
of each detected nucleus range from 2-4 kpc, which is about 3-6 times smaller than the average effective 
CO diameter of nearby spiral galaxies. This is strong evidence that the high concentration of molecular 
gas is the result of tidal dissipation in ongoing mergers. There is no correlation between the optical 
emission-line classification of the nuclei (i.e., Seyfert, LINER, or H II) and the presence of detectable 
molecular gas; however, there is a clear indication that the relative amount of molecular gas increases 
with the relative level of activity as measured via radio power and optical/near- infrared emission- line 
strength. Star formation rates are estimated to be in the range ~ 30 - 290 M year" 1 nucleus" 1 by 
making assumption that the radio and infrared emission arise from supernovae and dust heating by 
massive stars, respectively; corresponding gas consumption timescales are 1 — 7 x 10 7 years. The nuclei 
detected in CO are extremely red at near-infrared wavelengths, suggestive of much dustier environments 
than in the companions undetected in CO. Column density estimates are Ah 2 ~ io 24 ~ 25 cm~ 2 , which 
correspond to more than 1000 magnitudes of extinction toward the nuclei at visual wavelengths. Finally, 
the molecular gas mass densities and line-of-sight velocity dispersions show significant overlap with 
stellar densities and line-of-sight stellar velocity dispersions of local elliptical galaxies with My < — 19 
mag, including rapidly rotating ellipticals with disky isophotes and power-law light profiles as well as 
slowly rotating ellipticals with boxy isophotes and cores. This provides strong evidence that the CO-rich 
nuclei of these ULIGs have the phase-space density of gas necessary to form the stellar cores of elliptical 
galaxies. 

Subject headings: galaxies: active — galaxies: interacting — galaxies: ISM — galaxies: individual — 
ISM: molecules — infrared: galaxies - 



1. INTRODUCTION 

Recent surveys have provided tantalizing evidence that 
galaxy mergers at high redshift may be the precursors of 
local massive, quiescent galaxies. The fraction of both 
optically-selected mergers and submillimeter/infrared- se- 
lected galaxies (the latter of which are observed locally to 
be mergers: e.g. Joseph & Wright 1985) has been shown 
to increase with increasing redshift (Patton et al. 1997; 
Smail et al. 1997; Barger et al. 1998; Hughes et al. 1998; 



Le Fevre et al. 2000), and the space density of the high- 
redshift submillimcter population is observed to be similar 
to the space density of present day massive elliptical galax- 
ies (Trentham 2001). Further, N-body simulations have 
shown that the collision of disk galaxies naturally leads to 
the formation of massive elliptical galaxies (e.g., Barnes & 
Hernquist 1992 and references therein). 

Given the likely connection between galaxy mergers and 
massive elliptical galaxy formation, reconstruction of the 
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evolutionary steps - beginning with the initial interac- 
tion of the progenitors and ending with the final merger 
byproduct - will permit a determination of how the ac- 
tivity and dynamical state of galaxy mergers, and thus 
of massive galaxies, evolve as a function of evolutionary 
phase. Molecular gas is a key ingredient in the merger 
process; N-body simulations employing a dissipational gas 
component (i.e., Barnes & Hcrnquist 1996; Mihos & Hern- 
quist 1996) show that gas in the spiral disks of the progen- 
itor galaxies is gravitationally torqued, resulting in either 
stripping of gas from the disk, or a net inflow of material to 
the circumnuclcar regions of each progenitor. Multiwave- 
length observations of advanced infrared luminous mergers 
provide evidence that the enhanced molecular gas densi- 
ties lead to prodigious star formation and, in some cases, 
mass accretion onto supermassive nuclear black holes oc- 
curs (see review by Sanders & Mirabel 1996). 

The scenario presented above has been assembled from 
both theoretical and empirical work. However, observa- 
tional gaps exist, and feedback mechanisms (i.e., star for- 
mation, supernovae winds) are not understood well enough 
to implement in simulations. In order to fill in the ob- 
servational details regarding the role of molecular gas in 
the intermediate evolutionary stages of the most luminous 
mergers, a CO survey of ultraluminous infrared galaxies 
(ULIGs: L m [8 - IOOOmhi] £ 1O 12 L ) with double nu- 
clei and z < 0.15 has been initiated. These observations 
make use of the Owens Valley Radio Observatory (OVRO) 
Millimeter Array to resolve CO emission that may be as- 
sociated with each nucleus. This sample of ULIGs has 
been well-studied at other wavelengths (e.g., Condon et 
al. 1991, 1992; Murphy et al. 1996, 2001; Sanders, Scov- 
ille, & Soifer 1991; Sanders et al. 1988a, b; Scoville et al. 
2000; Soifer et al. 2000; Solomon et al. 1997; Surace 
et al. 1998, Surace & Sanders 1999; Surace, Sanders, & 
Evans 2000; Veilleux, Sanders, & Kim 1999); thus cor- 
relations between the dynamical state of the merger (as 
determined via optical-to-infrared imaging) and the stel- 
lar and AGN activity occuring (as determined via optical- 
to-infrared spectroscopy and radio interferometry) can be 
investigated. In addition, limits can be placed on the col- 
umn densities toward CO-detected nuclei and the gas mass 
density, which allow for a discussion of extinction in ULIGs 
and their possible association with stellar core formation in 
elliptical galaxies (e.g. Kormendy & Sanders 1992). Parts 
of our OVRO survey have been presented elsewhere (Evans 
et al. 1999; Evans, Surace, & Mazzarella 2000; Evans et 
al. 2000); galaxies discussed in these previous papers, as 
well as observations of related objects published by others 
(Sakamoto et al. 1999; Trung et al. 2001), are discussed 
along with new observations presented here for the first 
time. 

This paper is divided into six Sections. The selection cri- 
teria of the galaxy sample is discussed in §2. Section 3 is a 
summary of the new observations and data reduction. In 
§4, the CO(l — > 0) emission- line properties are presented, 
and the CO coordinates and morphologies are compared 
with the radio and HST optical and near-infrared mor- 
phologies. The method of calculating the molecular gas 
mass is briefly presented in §5. Section 6 contains a dis- 
cussion of the molecular properties of each nucleus in com- 
parison with the optical/near- infrared and radio emission 
lines, determinations of the star formation rate and the 



column density, and comparisons of the gas mass density 
to the stellar densities of elliptical galaxy cores. Section 
7 is a summary. Throughout this paper, Hq = 75 km 
s _1 Mpc _1 , qo = 0.5, and A = 0.0 are assumed. 

2. SAMPLE 

Four criteria for galaxies in the sample must be met in 
order to achieve the desired science goals. First, there must 
be sufficient information obtained at other wavelengths to 
show that they are strongly interacting merger pairs. Sec- 
ond, the nuclei must be far enough apart such that the 
progenitors are recognizable as separate galaxies. Third, 
the progenitors must also be far enough apart so that 
CO(l— >0) emission possibly present in the nuclei of both 
progenitors would be resolvable with OVRO (i.e, projected 
separations J> 2"). Fourth, the mergers must be well- 
studied at other wavelengths so that additional informa- 
tion (radio interferometry, optical spectroscopy, optical- 
to-infrared imaging) are available. 

Among the complete sample of 21 ULIGs from the 
IRAS Bright Galaxy Sample (BGS, with 60 flux den- 
sities, /60/jm > 5.24 Jy: Soifer et al. 1987) selected by 
Sanders et al. (1988a) and the IRAS Warm Galaxy Sam- 
ple (feofim > 1-5 Jy and .fc^m/iWm ^ 0.20: Sanders et 
al. 1988b), five galaxies fit the above criteria, and their 
general properties are summarized in Table 1. Three of 
these ULIGs (IRAS 08572+3915, IRAS 13451+1232, Mrk 
463) are classified as far-infrared "warm" galaxies, and two 
ULIGs (IRAS 12112+0305, IRAS 14348-1447) are classi- 
fied as far-infrared "cool" (i.e., /25 M m/ ' f&o^m < 0.20) galax- 
ies. 

3. OBSERVATIONS AND DATA REDUCTION 

Aperture synthesis maps of CO(l — * 0) emission from 
IRAS 08572+3915, Mrk 463, and IRAS 12112+0305 were 
made with the Owens Valley Radio Observatory (OVRO) 
Millimeter Array during five observing periods from 1996 
November to 1999 April (see Evans et al. 1999 and Evans, 
Surace, & Mazzarella 2000 for discussions of observations 
and data reduction of IRAS 13451+1232 and IRAS 14348- 
1447, respectively). The array consists of six 10.4m tele- 
scopes, and the longest observed baseline was 242m. Each 
baseline was configured with 120x4 MHz digital correla- 
tors. Observations done in low-resolution configuration 
provided a ~ 4'.'0 (FWHM) synthesized beam with natural 
weighting, and observations in high-resolution configura- 
tion provided a beam of ~ 2'/0 (FWHM). For the galax- 
ies IRAS 12112+0305 and Mrk 463, observations were 
also made with the equatorial configuration to improve 
the shape of the beam. During the observations, nearby 
quasars were observed every 25 minutes to monitor phase 
and gain variations, and 3C 345, 3C 454.3, and 3C 273 
were observed to determine the passband structure. Fi- 
nally, observations of Uranus were made for absolute flux 
calibration. The journal of observations is provided in Ta- 
ble 2. 

The OVRO data were reduced and calibrated using 
the standard Owens Valley array program MMA (Scov- 
ille et al. 1993). The data were then exported to the 
mapping program DIFMAP (Shepherd, Pearson, & Taylor 
1995) and to the National Radio Astronomy Observatory 
(NRAO) Astronomical Image Processing System (AIPS) 
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for extraction of spectra. 

4. RESULTS 

Radio data from the literature have been compiled for 
the ULIGs in the sample, and the radio source coordinates 
are listed in Table 3 along with the CO(l — ► 0) coordinates 
measured in this study. Both the CO and radio coor- 
dinates for each nucleus agree to within <; 0"5 accuracy. 
Hubble Space Telescope (HST) data of all five ULIGs have 
also been obtained from the HST archive. The optical and 
near-infrared morphologies of each galaxy are complex and 
will be discussed in the following subsections in relation to 
the CO distribution. The CO emission line properties are 
summarized in Table 4. 

4.1. IRAS 08572+3915 

Figure la shows a 0.8 /im image of the warm galaxy 
merger IRAS 08572+3915 taken with the HST WFPC2 
instrument (see Surace et al. 1998). The merging system 
consists of two spiral galaxies with a nuclear separation of 
5.6 kpc (5'.'4). A 10 kpc tidal tail extends northward of the 
northwestern galaxy (hereafter IRAS 08572+3915NW), 
and another tail extends for approximately 7 kpc eastward 
of the southeast galaxy (hereafter IRAS 08572+3915SE). 
Several unresolved knots arc visible along the southern 
tidal tail and within the inner 3 kpc of both nuclei. 

Figure lb shows the contours of the CO emission from 
IRAS 08572+3915 superimposed on a three-color near- 
infrared image of the merger obtain with the HST Near- 
Infrared Camera and Multi-Object Spectrometer (NIC- 
MOS) camera; 1.1, 1.6, and 2.2 (im are shown as blue, 
green, and red, respectively (see also Scoville et al. 2000). 
The NW nucleus is clearly much redder than the SE 
nucleus, and CO emission is detected only on the NW 
nucleus. Mid-infrared (12-25 jum) imaging of IRAS 
08572+3915 (Soifer et al. 2000) show all of the mid- 
infrared emission to be emanating from an unresolved com- 
ponent associated with the NW nucleus, indicating that 
the source of the strong mid-infrared luminosity is in the 
NW nucleus and is confined to a region less than 250 pc 
across. The CO emission is primarily confined to an unre- 
solved (2'.'1) component less than 2.2 kpc in extent, with a 
total CO flux of 11.8+1.4 Jy km s" 1 . This CO flux is con- 
sistent with the measurement of 9.0+1.8 Jy km s _1 made 
by Solomon et al. (1997) with the IRAM 30m telescope, 
and indicates that practically all of the CO emission mea- 
sured in the 22" beam of the 30m is associated with the 
NW galaxy. The difference between the single-dish and 
OVRO observations can be used to constrain the CO flux 
of the SE nucleus; the difference of —2.8+3.2 Jy km s _1 re- 
sults in an upper limit of 0.4 Jy km s _1 . However, this flux 
is 10 times less than the 3a upper limit of 4.2 Jy km s _1 
calculated using the root-mean-square (rms) noise in the 
OVRO data and assuming a velocity width equal to that 
of the NE nucleus. Therefore, the latter flux is adopted 
as the upper CO flux limit of IRAS 08572+3915SE (Table 
4). The CO(l — > 0) spectrum of NW nucleus is also shown 
in Figure lb. 

4.2. IRAS 12112+0305 

Figure 2a shows a 0.8^m HST archival image of the 
cool galaxy merger IRAS 12112+0305. Unlike IRAS 



08572+3915, the distinction between the progenitor galax- 
ies is less obvious at visible wavelengths. A tidal feature is 
observed to extend northward of the nuclear region, and a 
prominent tail extends to the west, then to the southeast 
of a compact knot in the southwest region of the nuclear 
complex. Various unresolved knots are observed through- 
out the entire system. 

Figure 2b shows the CO emission from IRAS 
12112+0305 superimposed on a three-color NICMOS im- 
age; 1.1, 1.6, and 2.2 /im are shown as blue, green, and red, 
respectively (see also Scoville et al. 2000). The complex 
nuclear structure of the WFPC2 image is reduced to two 
high surface brightness components separated by 3.7 kpc 
(2'.'9) in the NICMOS image; these components appear 
to be the nuclei of the progenitor galaxies. Unresolved 
(2.7 kpc = 2'.'2) CO emission is observed to be associated 
with both nuclei, with approximately 75% and 25% of the 
emission in the northeast (hereafter IRAS 12112+0305NE) 
and the southwest (hereafter IRAS 12112+0305SW) com- 
ponents, respectively. The total CO flux emitted from 
IRAS 12112+0305 is 39.7 Jy km s" 1 , within 15% of the 
CO measurement of 46 Jy km s _1 made with the NRAO 
12m (half-power beamwidth = 55": Sanders, Scoville, & 
Soifer 1991). A CO(l — > 0) spectrum of each nucleus is 
also shown in Figure 2b. 

4.3. IRAS 13451+1232 

The CO data for IRAS 13451+1232 (aka PKS 
1345+1232) are discussed in detail by Evans et al. (1999), 
and are therefore omitted here. 

4.4. Mrk463 

Figure 3a shows the 0.8^m HST archival image of the 
warm galaxy Mrk 463 (IRAS 13536+1836) from Surace et 
al. (1998). The two nuclei of the progenitor galaxies are 
3.7 kpc (4"1) apart, and both possess complex morpholo- 
gies which include multiple circumnuclear star clusters and 
faint tidal features. 

Figure 3b show the CO emission from Mrk 463 super- 
imposed on a three-color NICMOS image from Hines et 
al. (2002); 1.1, 1.6, and 2.07/xm are blue, green and red, 
respectively. The unresolved (2.2 kpc = 2"4) CO emis- 
sion is associated with the Seyfert 2 eastern nucleus (here- 
after Mrk 463E), which has been shown by Mazzarella et 
al. (1991) to have extremely red near-infrared colors con- 
sistent with a dust-enshrouded quasar (see also §6.5) as 
well as radio jet emission. Mrk 463E is also unresolved at 
2.07,um, indicating that the emission at this wavelength is 
confined to a region less than 200 pc in extent. The to- 
tal measured CO flux for Mrk 463E is 6.8+0.9 Jy km s"\ 
which agrees with the IRAM 30m measurement of 7.2+0.8 
Jy km s" 1 for Mrk 463 obtained by Alloin et al. (1992) 
to within 6%. The difference between the single-dish and 
OVRO observations can be used to constrain the CO flux 
of the W nucleus; the difference of 0.4 ± 1.7 Jy km s _1 re- 
sults in an upper limit of 2.1 Jy km s _1 . This is consistent 
with the 3a upper limit of 2.7 Jy km s _1 calculated using 
the rms noise in the OVRO data and assuming a velocity 
width equal to that of the E nucleus. The higher flux is 
adopted as the upper CO flux limit of Mrk 463W (Table 
4). The CO(l — > 0) spectrum of E nucleus is also shown 
in Figure 3b. 
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4.5. IRAS 14348-144 7 

The CO data for IRAS 14348-1447 are discussed in de- 
tail by Evans, Surace, & Mazzarella (2000), and thus the 
discussion will be omitted here. 

5. CO LINE LUMINOSITIES & MOLECULAR GAS MASSES 

Table 4 summarizes the emission-line properties of the 
five ULIGs. For a A = Universe, the luminosity distance 
for a source at a given redshift is 

D L = cH^q 2 {zq Q + (qo - 1) (y/^Qoz + 1 - l) } [Mpc]. 

(1) 

Given the measured CO flux, SVjoAu [Jy km s 1 ], the CO 
luminosity of a source at redshift z is 



CO 



2^o 2 bs 



ScoA«D£(l + z)- 3 [K km s- 1 pc 2 ], 



(2) 

(Solomon, Downes, & Radford 1992) where c [km s x ] is 
the speed of light, k [J K" 1 ] is the Boltzmann constant, 
and ^bs [Hz] is the observed frequency. In terms of useful 
units, £co(i->o) can ^ e wr ^ ten as 



L' co = 2.4 x 10 3 



Jy km s~ 



[K km s" 1 pc 2 ] 



(3) 

To calculate the mass of molecular gas in these galax- 
ies, the reasonable assumption is commonly made that the 
CO emission is optically thick and thermalizcd, and that 
it originates in gravitationally bound molecular clouds. 
Thus, the ratio of the H 2 mass and the CO luminosity 
is given by 



M(H 2 ) 
a = — — <x 



L 



CO 



[M^Kkms- 1 pc 2 )- 1 ], (4) 



where n(H 2 ) and Tb are the density of H2 and bright- 
ness temperature for the CO(l— >0) transition (Scoville & 
Sanders 1987; Solomon, Downes, & Radford 1992). Mul- 
titransition CO surveys of molecular clouds in the Milky 
Way (e.g. Sanders et al. 1993), and in nearby starburst 
galaxies (e.g. Glisten et al. 1993) have shown that hotter 
clouds tend to be denser such that the density and temper- 
ature dependencies cancel each other. The variation in the 
value of a is approximately a factor of 2 for a wide range 
of kinetic temperatures, gas densities, and CO abundance 
(e.g. a = 2 — 5M Q [K km s -1 pc 2 ] -1 : Radford, Solomon, 
& Downes 1991). Solomon et al. (1997) and Downes & 
Solomon (1998) have made use of dynamical mass esti- 
mates of a low-redshift infrared galaxy sample observed 
in CO with the Plateau de Bure Interferometer to argue 
that a may, in some cases, be as low as 1 M Q (K km s -1 
pc 2 )" 1 . 

In order to determine an appropriate value of a for the 
present sample of galaxies, the assumption will be made 

7 The assumption has been made that both nuclei of a given ULIG 



that the mass of the inner region of each galaxy is primar- 
ily comprised of molecular gas. Thus, the molecular gas 
mass will be assumed to be equal to the dynamical mass, 



M(H 2 ) ~ M dy 



A^FWHM^CO 



G 



226 



/ AfFWHM \ ' 

\ km s -1 / 



[M 



(5) 



where Avfwhm is the full CO velocity width at half the 
maximum flux density (Table 4) and i?co is the radius of 
the CO distribution. Given that the CO emission associ- 
ated with each nucleus was not resolved with the OVRO 
observations (Figure 1-3), i?co is estimated in two ways. 
First, the assumption is made that the molecular gas is 
similar in extent to the radio emission (Condon et al. 1991: 
Table 5); this is plausible given that the radio emission 
may, with the exception of IRAS 13451+1232 and Mrk 
463 (sec §6.2), be due to supernovae from star formation 
that has been fueled by the molecular gas (hereafter, size 
estimate 1). Second, the CO extent is estimated by as- 
suming optically thick, thermalized gas with a unity filling 
factor and a blackbody temperature, Tbb, equal to that of 
the dust 7 (hereafter, size estimate 2). Given the latter 
assumption, Tbb can be determined using the 60^m and 
lOO^m flux densities 



tbb 



-(1 + *) 



82 



m (0.3/60 ( um//l00 ( um) 



0.5 



[K] (6) 



(see discussion in Solomon et al. 1997). Thus, the radius 
of the CO distribution is calculated via 



i?co = 



7 CO 



TrTbbAvpwHM 



[PC] 



(7) 



(see Table 5). The sizes estimated via these two methods 
are with a factor of 2-4 of each other, with the size esti- 
mate 1 yielding smaller CO extents. Adopting the more 
conservative CO extent determined from size estimate 2 
yields M dyn /L' co = 1.0 - 2.5 M Q (K km s" 1 pc 2 )" 1 with 
a mean value of 1.6 M (K km s" 1 pc 2 )" 1 . These val- 
ues are similar to Mdyn/L'co = 0-8 — 2.5 M Q (K km s _1 
pc 2 ) -1 derived by Solomon et al. (1997) for a sample of 
37 ULIGs; the average of their sample is 1.3 M Q (K km 
s" 1 pc 2 )" 1 . We adopt a of 1.5 M Q (K km s" 1 pc 2 )~\ 
which is lower by a factor of 3 than the value determined 
for the bulk of the molecular gas in the disk of the Milky 
Way (Scoville & Sanders 1987; Strong et al. 1988). The 
derived molecular gas masses are listed in Table 4. 

6. DISCUSSION 

The data presented in Figures 1-3 clearly illustrate the 
advantage of these new observations. The beam size is a 
factor of 10-30 smaller than previous single-dish CO ob- 
servations of these double-nucleus galaxies (Sanders, Scov- 
ille, & Soifer 1991; Alloin et al. 1992; Sanders et al. 1993; 
Solomon et al. 1997), and the improved resolution is suf- 
ficient to determine the location of the CO emission rel- 
ative to both the stellar distribution and radio emission 

have the same dust temperature, which may be incorrect. 
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(§4). The beam size is also sufficient to show that the up- 
per limit of the CO extent of each nucleus (2.1-4.4 kpc) 
is significantly less than the average effective CO diameter 
of 11.5 ± 7.5 kpc determined for a sample of ~ 140 nearby 
spiral galaxies (Young et al. 1995; see also Evans, Surace, 
& Mazzarella 2000). 8 This latter comparison may be an 
indication that the molecular gas in the progenitor galax- 
ies, most of which appear to be spiral galaxies based on 
the presence of spiral and/or tidal features (Figures 1-3; 
see also Figure 1 of Evans, Surace, & Mazzarella 2000), 
has already been driven inwards via the gravitational in- 
teraction of the galaxy pairs. 

The most obvious characteristic of these data is that 
not all of the nuclei of these intermediate-stage ULIGs are 
detected in CO(l — > 0). This issue has been previously 
discussed; Evans et al. (1999) presented CO observations 
of the warm ULIG IRAS 13451+1232 which show molecu- 
lar gas to be only detected in radio-loud AGN nucleus, and 
previous CO observations of the cool ULIGs Arp 220 and 
IRAS 14348-1447 show molecular gas associated with both 
nuclei (Sakamoto et al. 1999; Evans, Surace, & Mazzarella 
2000). The new data presented here for the additional 
galaxies support that general trend; both of the warm 
ULIGs (IRAS 08572+3915 and Mrk 463) have molecular 
gas detected only on the redder, radio luminous nucleus, 
whereas molecular gas is detected on both nuclei of the 
cool ULIG IRAS 12112+0305. This effect is not entirely 
due to differences in the sensitivities of the observations; 
two of the three warm ULIGs (IRAS 08572+3915 and Mrk 
463) are the closest objects in the sample, and their un- 
detected nuclei have L' co upper limits that are a factor of 
3.6-7.5 less than the CO luminosities of the least CO lumi- 
nous cool galaxy nucleus (IRAS 12112+0305SW). In terms 
of molecular gas mass, IRAS 08572+3915SE and Mrk 
463W have molecular gas masses a factor of 3.0-6.2 less 
than that of the Milky Way Galaxy (M MW [H 2 ] ~ 2.5 x 10 9 
M©). Note, however, that while the single-dish and OVRO 
observations yield comparable flux values, consistent with 
the possibility that the molecular gas is associated only 
with the more radio luminous nucleus, it cannot be ruled 
out in the cases of IRAS 13451+1232 and Mrk 463 that 
their galaxy pairs have L' co ratios similar to that of the 
cool ultraluminous galaxy pairs (see Table 6). 

The wealth of multiwavelength data available for this 
sample make it possible to explore correlations between 
the molecular gas and both the activity occuring in each 
nucleus and the estimated stellar masses. In addition, sim- 
ple assumptions about the true extent of the molecular gas 
will be made to derive crude determinations of the column 
density and gas mass density of the detected nuclei. These 
issues, and their relevance to ULIGs, will be explored for 
the majority of the discussion section. 

6.1. Emission- Line Classification and Line Strength 

Table 1 contains a compilation of the optical emission 
line classifications of each nucleus (where available) . Such 
diagnostics have been used over the past two decades to 
classify galaxies as Seyfert galaxies, Low Ionization and 
Nuclear Emission Region galaxies (LINERs), and HII- 
region like galaxies. The Seyfert and HII region-like 



classification are indicative of AGN and starburst activ- 
ity, respectively, whereas the LINER classification may 
be due to shocks produced by either an AGN with a 
shallow power law, or by supernovae events. Extensive 
studies of ULIGs and Luminous Infrared Galaxies (LIGs: 
L IR = iQii-OO-n.99 Lq ) ghow that the f raction f HII 

region-like galaxies decreases with increasing infrared lu- 
minosity, the fraction of Seyfert galaxies increases, and 
the fraction of LINERs stays roughly constant (e.g., Kim 
& Sanders 1998). 

Both Seyfert nuclei (classified as Type 2 optically, but 
Type 1 via polarized or infrared emission: e.g., Veilleux, 
Kim & Sanders 1997) contain copious amounts of molecu- 
lar gas, whereas only 4 of the 6 LINER nuclei have detected 
CO. This is in contrast to what has been observed in dy- 
namically young, widely separated (i.e., > 20 kpc) ULIGs 
(Trung et al. 2001) in which all Seyfert and LINER nuclei 
in their sample have 0.5 — 1.6 x 10 10 M Q of molecular gas 9 
associated with them. If indeed the galaxies identified by 
Trung et al. are components of ultraluminous mergers, 
then either there is a difference between ULIGs comprised 
of widely separated galaxy pairs and intermediate stage 
ULIGs (e.g., gas may be stripped from one of the galaxies 
after the initial encounter and transferred to its compan- 
ion as the merger progresses, thus leaving some relatively 
low-luminosity LINER nuclei gas-poor), or the samples are 
statistically too small to derive firm conclusions. 

Another obvious avenue to explore is a possible correla- 
tion between the amount of molecular gas and the strength 
of optical and near-infrared emission lines. In other words, 
it is reasonable to expect that emission lines associated 
with the gas-rich nuclei are more luminous because there 
is more gas present to fuel more vigorous activity. Recom- 
bination lines of hydrogen such as Paschen a and Ha are 
suitable for such a comparison; these photons are a direct 
probe of ionizing photons from young, massive stars and 
AGN. 

Paschen a measurements exist for both nuclei of IRAS 
08572+3915, IRAS 12112+0305, and IRAS 14348-1447 
(for the others, only Paa data of the radio-loud AGN nu- 
clei are published) and the ratio of Paa flux density of 
nuclei in the ULIGs sample are tabulated along with CO 
flux density ratios in Table 6. With the exception of IRAS 
12112+0305, the data appear to confirm that, at 1.9^m, 
more luminous line emission is associated with more molec- 
ular gas. The Ha measurements (Table 6), which suffer 
~ 10 magnitudes more extinction than Paa and are thus 
less reliable indicators of relative nuclear activity, are also 
consistent with this result. 

6.2. Radio Data 

Another tracer of activity in galaxies is radio continuum 
emission. In galaxies, it is primarily due to synchrotron 
emission from AGN or supernovae. In the latter case, the 
emission is an indication of recent high-mass star forma- 
tion (i.e., stars with masses 8 M Q , which live <; 3 x 10 7 
years: Condon 1992). 

Two of the galaxy pairs (Mrk 463 and IRAS 
13451+1232) in the sample clearly have radio emission due 



8 Note that the CO luminosities measured with OVRO arc consistent with prior singlc-dish measurements (see §4), indicating that there 
cannot be significant extended low surface brightness CO emission. 

9 The molecular gas masses in Trung et al. (2001) have been recalculated with a = 1.5 Mq (K km s _1 pc 2 ) -1 . 
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to an AGN; the emission is in the form of radio jets ema- 
nating from the nuclear region (e.g. Mazzarella et al. 1991; 
Stanghcllini et al. 1997). The remaining three galaxy pairs 
have fairly compact radio emission (Condon et al. 1990, 
1991), with upper limit sizes on the order of 100-200 pc. 

The ULIGs IRAS 14348-1447 and IRAS 12112+0305 are 
the only two ULIGs for which both nuclei are detected 
in CO and for which no obvious evidence of radio emis- 
sion due to AGN presently exists. As is the case with the 
optical/near-infrarcd line emission, more luminous radio 
emission is associated with more molecular gas. 

6.3. Stellar Masses 

It is commonly asserted that ultraluminous infrared 
galaxies are created by the merger of equal-mass disk 
galaxies. In order to estimate the relative masses of galax- 
ies in each ultraluminous galaxy pair so that a compari- 
son can be made with their relative molecular gas masses, 
the relative _B-band, /-band and //-band luminosities of 
the component galaxies in the five ULIGs have been mea- 
sured within polygonal regions that encompass the emis- 
sion from each galaxy on HST and ground-based images 
published by Surace et al. (1998), Surace, Sanders, & 
Evans (2000), and Scoville et al. (2000). The //-band flux 
density ratios among the paired nuclei of all five ULIGs 
are in the range 0.8-2.5, with Mrk 463E-to-W having the 
highest ratio. Since the dust-enshrouded quasar in Mrk 
463E contaminates the //-band light much more so than 
the B-band light (Mazzarella et al. 1991), the stellar mass 
ratio traced by the //-band flux ratio of Mrk 463E to Mrk 
463W is an upper limit. IRAS 08572+3915SE and IRAS 
13451+1232NW also have very luminous, red nuclei con- 
sistent with dust-enshrouded QSOs (see §6.5). Therefore, 
the /-band and _B-band flux ratios are more relevant than 
the //-band ratios for comparing the merging galaxies with 
minimal contamination from their nuclei. The /-band and 
£>-band flux density ratios among the paired nuclei of all 
five ULIGs are within the range 0.6-1.4, with the lumi- 
nosities of IRAS 08572+3915SE, IRAS 13451+1232NW, 
and Mrk 463W being, within the uncertainties, essentially 
the same as their CO-rich companions. Comparable stel- 
lar masses are also implied by the similar disk sizes of the 
individual merging galaxies (e.g., see Figures 1-3). 

The range in optical flux density ratios compared to the 
ratios of molecular gas masses for galaxies in each pair (Ta- 
ble 6) implies that there is not a constant ratio of L' co /Lb- 
This is consistent with observations of a volume-limited 
sample of nearby spiral galaxies (Sage 1993); an examina- 
tion of L' co ratios for spirals in that sample paired such 
that their Lq ratios are in the range observed in these 
ULIRGs (0.6 - 1.4) shows notable scatter. Indeed, the 
Sage data show that the ULIG nuclei undetected in CO 
may be just below the threshold of the current observa- 
tions. 

6.4. Star Formation Rates and Gas Consumption 
Timescales 

Ultraluminous infrared galaxies provide one of the best 
samples of galaxies with which to study star formation in 
extreme environments. Several different diagnostics exist 
for calculating the star formation rate (SFR) in galaxies, 
and four diagnostics will be applied to the individual nuclei 



of the present sample. The calculated SFRs are tabulated 
in Table 7. 

6.4.1. SFR from Far-Infrared Luminosity 

The far-infrared luminosity is not an ideal technique for 
measuring the SFR of the present sample of galaxies due 
to (i) the ambiguity in what fraction of the far-infrared 
luminosity emanates from each nucleus and (ii) the un- 
known fraction of dust heating by AGN. Regardless, this 
technique will provide a point of comparison for the other 
diagnostics. The assumptions associated with this tech- 
nique are that (i) the infrared luminosity is primarily due 
to dust heated by an imbedded starburst, (ii) the optical 
depth to the starburst is high enough such that the in- 
frared luminosity is equal to the total luminosity of the 
starburst, (Hi) the starburst is continuous, and (iv) that 
the initial mass function (IMF) is well approximated by a 
Salpeter IMF. The SFR equation is 



SFR= 1.76 x 10 



-10 



L 



FIR 



\Mc. 



(8) 



(e.g. Scoville & Young 1983; Gallagher & Hunter 1986; 
Kcnnicutt 1998a), where Lfir (=L[40 — 500/zm]; see Table 
1) is calculated using the 60^m and 100/im fluxes densities 
and thus avoids the contribution to the infrared luminosity 
due to hot dust that may be heated by an AGN compo- 
nent (see Sanders & Mirabel 1996). The calculated SFRs 
are in the range 30-300 Mq yr _1 (see Table 7); however, 
the actual star formation rates may be significantly lower 
if AGN contribute significantly to the heating of dust in 
these ULIGs. 

6.4.2. SFR from Radio Luminosity (1.4 GHz) 

Condon et al. (1990, 1991) present 1.49 GHz and 8.4 
GHz radio flux densities and sizes of galaxies in the IRAS 
Bright Galaxy Sample. To calculate the SFR based on 
the 1.4 GHz luminosity, (i) a Salpeter IMF is assumed, 
(ii) the thermal Bremsstrahlung contribution to the 1.4 
GHz emission is assumed to be negligible, and (Hi) the 
assumption is made that AGN do not contribute signifi- 
cantly to the radio emission. Given the latter assumption, 
IRAS 13451+1232 and Mrk 463, both of which have radio 
jets associated with AGN, have been omitted from this 
calculation. The SFR equation is 



SFR= 1.36 x 10~ 28 



L\A GHz 



ergs s" 



Hz- 



[Moyr- 1 ] (9) 



(see Condon 1992; Haarsma et al 2000; Sullivan et al. 
2001). The resultant SFRs for IRAS 08572+3915NW 
and the nuclei of IRAS 12112+0305 and IRAS 14348- 
1447 are in the range 60-350 M Q yr" 1 . The SFR ratios 
SFR(1.4 GHz) / SFR(Lfir) are ~ 0.3, 1.2, and 1.8 for IRAS 
08572+3915, IRAS 12112+0305, and IRAS 14348-1447, 
respectively. The far-infrared and radio data thus give 
consistent estimates, within a factor of 3, for the SFR of 
each ULIG. 

6.4.3. SFR from Molecular Gas Surface Density and 

Extent 

The global Schmidt law for galaxies is 
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S SFR = (2.5 ± 0.7) x 10" 



( S g 



Vm qP c- 



1.4±0.15 



[M Q yr 1 kpC 



(10) 



(Kcnnicutt 1998a, b), where S gas is the molecular gas sur- 
face density (= M (H2) / V i? co ) , Ssfr is the SFR surface 



density, and a = 4.2 M c . 



value of a to 1.5 M (K km s 1 pc 2 
equation to be in terms of M Q yr -1 gives 



(K km s 1 pc 2 ) 1 . Changing the 
1 and modifying the 



SFR= (2.7 ±0.2) x 10" 



M pc- 



1.4±0.15 



CO 



V kpc 



[M yr- 1 ], (11) 

where Dqo is the extent of the CO emission. Using the 
CO size estimates 1 and 2 discussed in §5 to determine 
E gas , SFR of 350-2110 M Q yr" 1 nucleus" 1 and 60-1080 
M Q yr -1 nucleus -1 , respectively, arc derived. These SFRs 
are factors of 4-8 and 2-3, respectively, higher than those 
derived from the 1.4 GHz luminosity. The SFRs are well 
within the uncertainties of equation 11; for example, low- 
ering the exponential term from 1.4 to 1.25 (size estimate 
1) or 1.3 (size estimate 2) yields SFR consistent with the 
radio luminosity determinations of SFR. 

6.4.4. SFR from Molecular Gas Mass and Gas 
Consumption Timescale 

Using a sample of 36 infrared-selected starburst galax- 
ies, Kennicutt (1998b) calculated the ratio of their molec- 
ular gas masses to their star formation rates derived from 
Lfir to derive gas consumption timescales. The average 
gas consumption timescale is tdyn ~ 4.6 x 10 7 years; us- 
ing a = 1.5 Mq (K km s" 1 pc 2 ) -1 instead of 4.2 M Q 
(K km s" 1 pc 2 )" 1 yields a timescale of idyn ~ 1-6 x 10 7 
years; this is significantly less than idyn for normal spiral 
galaxies, which are in the range 1 — 8 x 10 8 years (Ken- 
nicutt 1998b). The shorter gas consumption timescale of 
infrared starburst galaxies is further confirmed by deriving 
£dyn for a sample of 37 ULIGs observed with the IRAM 
30m telescope by Solomon et al. (1997); using their tab- 
ulated Lfir/L' co ratios and ana = 1.5 M (K km s" 1 
pc 2 ) -1 , the average ta yn is 4.7 x 10 7 years. 

Table 7 lists the gas consumption timescales derived for 
the present sample of galaxies by making use of the SFR 
derived from both Lfir and L1.4 gh z - The average £dyn 
are 4.1 x 10 7 and 3.7 x 10 7 years, respectively, within 15- 
22% of tdyn derived for the Solomon et al. (1997) sample. 
By adopting td yn ~4x 10 7 years, the SFR based on tdyn 
and M(H 2 ) can be written as 



SFR = 2.5 x lO" 8 ( 



4 x 10 7 yrs 



Aid 



[M yr" 



(12) 

The SFRs derived from equation 12 are listed in Table 7. 

To summarize, the various techniques have been used 
to determine that the star formation rate for most of the 
CO-luminous nuclei is >100 M yr" 1 . As is expected for 
such dynamically evolving systems, these rates are consid- 
erably higher than the passive star formation rate of ~ few 
M yr -1 for the Milky Way galaxy. 



6.5. Extinction and Column Densities 

While optical and mid-infrared emission line diagnos- 
tics have been used to characterize the energy sources in 
ULIGs, it has long been known that significant amounts of 
extinction render such analysis inconclusive for imbedded 
nuclear sources. The effects of extinction in the molec- 
ular gas-rich nuclei of ULIGs can readily be seen by a 
simple analysis of near-infrared photometry. Color infor- 
mation in the near-infrared wavelength range are used to 
constrain the stellar population age and to provide some 
insight on how dusty the galaxies may be. Figure 4 is a 
plot of the difference in the 1.1 /im and 1.6 pja magni- 
tudes, mi. 1-1. 61 versus the difference in the 1.6 /xm and 
2.2 /im magnitudes, mi. 6 -2.2, for nuclei of ULIGs in the 
sample; data for the gas-rich nuclei of Arp 220 are also 
plotted for comparison. These data are plotted along with 
both an instantaneous burst Salpeter IMF which traces 
the evolution of the near-infrared colors of stars in the 
mass range 0.1 to 125 M up to an age of 0.3 Gyr, and 
for the near-infrared colors of QSOs with various contribu- 
tions of hot dust. (See Scoville et al. 2000 for an extensive 
discussion on NICMOS infrared color diagrams.) As is 
clear from this figure, (i) dust contributes to various de- 
grees to the near-infrared colors of all of the nuclei and 
(ii) the nuclei which have been detected in CO are signifi- 
cantly redder than the nuclei that have not been detected 
in CO. This provides reassuring proof that "red" nuclei are 
red at near-infrared wavelengths because of the presence 
of gas and dust, and also is an indication that light from 
the energy sources responsible for the emission observed 
at near-infrared wavelengths in these red nuclei is heavily 
affected by dust. 

An estimate of the optical extinction can be derived 
from the gas column densities. Column density limits have 
been derived from the molecular gas mass and by using 
both size estimates (see Table 8). These estimates trans- 
late into visual extinctions of Ay > 1000 magnitudes. 



6.6. Molecular Gas Mass Densities and Elliptical 
Galaxies Stellar Cores 

Kormendy & Sanders (1992) compared the stellar and 
molecular gas mass densities and velocity dispersions of 
a few ULIGs with the stellar mass densities and velocity 
dispersions of elliptical galaxy cores in an attempt to es- 
tablish an evolutionary connection between the two galaxy 
types. As a follow-up to their work, molecular gas mass 
densities have been calculated for each nucleus using both 
CO size estimations, and the results are compiled in Table 
8. The latter two size determinations give mass densities 
primarily in the range p ~ 10 2-4 M pc" 3 and are plotted 
along with the line-of-sight CO velocity dispersions, a, in 
Figure 5 (i.e., the cooling diagram). Given the obvious 
uncertainties in determining p and a, the area of Figure 
5 occupied by the ULIGs should be treated as the gen- 
eral range of density and dispersion values that interme- 
diate stage mergers may possess. The stellar mass density 
and line-of-sight velocity dispersions of nearby elliptical 
galaxies, obtained from Faber et al. (1997), are also plot- 
ted. The gas mass densities and velocity dispersions of the 
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ULIGs are consistent with the stellar mass densities and 
dispersion of elliptical galaxies with absolute visual mag- 
nitudes My < —19. The ULIG data points overlap the 
regions occupied by both rapidly rotating elliptical galax- 
ies with disky isophotcs and power-law light profiles, which 
Faber et al. (1997) speculate may be the end-products of 
gaseous merger events, and slow rotating elliptical galaxies 
with boxy isophotes and cores. From the data plotted in 
Figure 5, it can be concluded that the phase-space density 
of star-forming molecular gas is sufficient to form the stars 
in the nuclear regions of luminous elliptical galaxies. 

7. SUMMARY & CONCLUDING REMARKS 

Molecular gas observations of a sample of double- 
nucleus ultraluminous infrared galaxies (ULIGs) were pre- 
sented. This observations were motivated by the dearth 
of CO(l — > 0) data for ULIGs at intermediate stages of 
evolution (i.e., with nuclear separations of 3-5 kpc). The 
following conclusions have been reached: 

(1) Three observed double-nucleus objects with warm far- 
infrared dust temperatures have detected CO emission as- 
sociated only with one of the two nuclei - the nucleus that 
is reddest and has the highest radio continuum emission. 
Two observed objects with cool far-infrared dust temper- 
atures have detected CO emission associated with both 
nuclei. Molecular gas masses for the detected nuclei are 
estimated to be in the range 0.1 - 1.2 x 10 10 M Q . 

(2) The relative amount of molecular gas in each galaxy 
pair appears to be correlated with the relative levels of 
activity as measured by both optical/near- infrared recom- 
bination line emission and radio flux density. The presence 
of LINER and Seyfert nuclei combined with the high cen- 
tral concentration of CO provides evidence that molecular 
gas is an important component in fueling AGN and vigor- 
ous, massive star formation. 

(3) Where possible, star formation rates have been com- 
puted for each nucleus using four different techniques. 
Star formation rates estimated to be 20-290 M Q yr _1 
nucleus -1 , with corresponding gas consumption timescales 
of 1 — 7 x 10 7 years. The star formation rates may be 
considerably lower (and the gas consumption timescales 
considerable longer) if AGN contribute significantly to the 
infrared and radio emission observed in these galaxy pairs. 

(4) The nuclei with associated molecular gas appear to 
have significantly redder infrared colors than their com- 
panions that lack CO detections, suggestive of dustier en- 
vironments in the CO-luminous nuclei. 

(5) Column densities for the nuclei of the ULIGs were esti- 
mated to be ~ 10 24 ~ 25 cm~ 2 , which corresponds to greater 
than 1000 magnitudes of visual extinction. 

(6) Where possible, molecular gas mass densities were es- 
timated for each nucleus. A reasonable estimate for the 
range of molecular gas densities is 10 2-4 M Q pc~ 3 . Such 
values, combined with measured line-of-sight CO(l — > 0) 
velocity dispersions, are equivalent to the stellar mass den- 
sities and velocity dispersions of elliptical galaxies with 
My < — 19 and indicate that the molecular gas has a suf- 
ficiently high phase-space density to form the stars in el- 
liptical galaxy cores. 

The present sample is small (5 galaxies), but the analysis 
presented here shows compelling results that may be solid- 
ified or refuted with a larger, more statistically significant 
sample of ULIGs in their intermediate phases of evolution. 



A CO(l — ► 0) survey is presently underway to observe all 
of the intermediate stage ULIGs in the IRAS 2 Jy sample. 
These data will be combined with multiwavelength data in 
a manner similar to what has been presented for the five 
ULIGs in this paper. Further improvement of the present 
dataset will also be possible with high-resolution CO ob- 
servations with the soon-to-be commissioned Smithsonian 
SubMillimeter Array (SMA) and the upcoming Combined 
Array for Research in Millimeter Astronomy (CARMA); 
these interferometers will make it possible to obtain CO 
resolutions approaching that of the 2[im HST data, thus 
removing the ambiguity in the true extent of the molecular 
gas. 
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CO IN DOUBLE NUCLEI ULIGS 



Figure Captions 



Figure 1. a) Hubble Space Telescope WFPC2 image of IRAS 08572+3915. b) CO(l -> 0) contours of the merger 
superimposed on a three-color composite NICMOS image (blue=l.l /im, green = 1.6 /im, red=2.2 /im). The CO contours 
are plotted as 30%, 40%, 50%, 60%, 70%, 80%, 90%, and 99% the peak flux of 8.4 Jy km s" 1 beam" 1 . The CO emission 
is unresolved, with a beam FWHM of 2'.'4 x l'/8 at a position angle of 77.9°. The apparent northwest extension of CO 
emission may be an artifact. The CO(l— >0) spectrum of the NW nucleus is also shown. For both images, north is up and 
east is to the left. 

Figure 2. a) Hubble Space Telescope WFPC2 image of IRAS 12112+0305. b) CO(l -> 0) contours of the merger 
superimposed on a three-color composite NICMOS image (blue=l.l /im, green = 1.6 /im, rcd=2.2 /im). The CO contours 
are plotted as 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, and 99% the peak flux of 24.7 Jy km s' 1 beam" 1 . The CO 
emission for each nucleus is unresolved, with a beam FWHM of 2? 7 x l'/9 at a position angle of -66.8°. The CO(l— > 0) 
spectra of the SW and NE progenitors are also shown. For both images, north is up and east is to the left. 
Figure 3. a) Hubble Space Telescope WFPC2 image of Mrk 463. b) CO(l — ► 0) contours of the merger superimposed on 
a three-color composite NICMOS image (bluc=l.l /im, green = 1.6 /im, rcd=2.07 /im). The 2/im emission is unresolved 
for Mrk 463E (the speckle pattern is an artifact of the NICMOS PSF), while Mrk 463W is clearly resolved. The CO 
contours are plotted as 50%, 60%, 70%, 80%, 90%, and 99% the peak flux of 6.4 Jy km s" 1 beam" 1 . The CO emission is 
unresolved, with a beam FWHM of 2'.'7 x 2'.'2 at a position angle of -61.9°. The CO(l— >0) spectrum of the E nucleus is 
also shown. For both images, north is up and east is to the left. 

Figure 4. mi. 1-1.6 vs. mi. 6-2. 2 color-color diagram for the five ULIGs; data for the nuclei of Arp 220 are also plotted. 
Fluxes are measured in a l'.'l diameter aperture centered on each nucleus. The locus for the evolution of an instantaneous 
starburst with a Salpeter IMF is also shown. The typical colors of optical PG QSOs are shown plotted with variable 
percentages of 2.2 /im emission due to warm dust. Lastly, a reddening vector based on the extinction curve of Rieke & 
Lebofsky (1985) modified at the shortest wavelengths is shown for two cases - a foreground screen of dust and a model in 
which the dust is uniformly mixed with the emitting sources. Adapted from Scoville et al (2000). 

Figure 5. A plot of the gaseous mass density versus the line-of-sight CO velocity dispersion for the five ULIGs listed 
in Table 1. The CO velocity dispersion, a, has been estimated by treating each CO line profile as Gaussian. Thus, 
a = Avfwhm/2.354. The small and large circles connected by lines represent mass density ranges limited using CO size 
estimates 1 and 2, respectively (see §5). Also plotted are the stellar mass densities and line-of-sight velocity dispersions 
of elliptical galaxies with both power-law light profiles and cores. The latter data are taken from Faber et al. (1997). 



Table 1 
The Sample 



Source 






Optical 
Class 


log L Fm 
log (L©) 


log Im 
log (L ) 


log Lboi 
log (L s ) 


/25//6O 


Nuclear 
Separation 


IRAS 08572+3915 
IRAS 08572+3915NW 
IRAS 08572+3915SE 


1021 


0.05821 
0.05821 


LINER 
LINER 


12.00 


12.10 


12.18 


0.23 


5" 4 


IRAS 12112+0305 
IRAS 12112+0305NE 
IRAS 12112+0305SW 


1246 


0.07268 
0.07268 


LINER 
LINER 


12.19 


12.26 


12.29 


0.06 


2"9 


IRAS 13451 + 1232 
IRAS 13451 + 1232NW 
IRAS 13451 + 1232SE 


1935 


0.1224 


Scyfcrt 2 


12.00 


12.25 


12.41 


0.35 


2"0 


Mrk 463 
Mrk 463E 
Mrk 463W 


905 


0.05087 
0.05087 


Seyfert2 

Seyfert2/LINER 


11.27 


11.76 


12.02 


0.72 


4'.'1 


IRAS 14348-1447 
IRAS 14348-1447NE 
IRAS 14348-1447SW 


1385 


0.0823 
0.0827 


LINER 
LINER 


12.23 


12.27 


12.30 


0.07 


3"4 



References. — Rcdshifts: Kim & Sanders 1998, except Mrk 463 (Mazzarclla et al. 1991). Emission Line Classification: 
Vcillcux, Kim, & Sanders 1999, except Mrk 463 (Chatzichristou & Vandcrricst 1995). Infrared Luminosities: Kim Sz Sanders 
1998, except Mrk 463 (Lir: Sanders et al. 1989; Lfir: using the equation for calculating Lfir given the luminosity distance 
and the IRAS flux densities provided in Sanders & Mirabel 1996). Infrared flux densities: Kim & Sanders 1998, except Mrk 
463 (IRAS Faint Source Catalogue). 
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Table 2 
Journal of Observations 



Phase Calibrator 



Source 


Dates 


JName 


Phase Coord. J2000.0 


3mm Flux Density 


Passband 








RA 


Dec 


(Jy) 


Calibrator 


IRAS 08572+3915 


1999 March 23 


.10927+390 


09:27:03.00 


39:02:20.59 


4.5 


3C 273 




1999 March 27 














1999 April 17 












IRAS 12112+0305 


1999 April 1 


3C 273 


12:29:06.66 


02:03:08.41 


17 


3C 273 




1999 April 12 












Mrk 463 


1996 November 9 


J1415+133 


14:15:58.78 


13:20:23.61 


1.5 


3C 273 




1996 December 25 










3C 345 




1997 January 12 










3C 454.3 



Table 3 

Measured Radio and CO Positions for the Detected Nuclei 



Source 


Radio coords .12000.0 
RA Dec 


CO coords J2000.0 
RA Dec 


ARA 

(") 


ADcc 
(") 


IRAS 08572+3915NW 


09:00:25.37 


39:03:54.16 


09:00:25.38 


39:03:53.82 


-0.12 


0.34 


IRAS 12112+0305NE 
IRAS 12112+0305SW 


12:13:46.02 
12:13:45.91 


02:48:41.33 
02:48:38.93 


12:13:46.03 
12:13:45.95 


02:48:41.33 
02:48:39.20 


-0.15 
-0.60 


0.00 
-0.27 


IRAS 13451 + 1232NW 


13:47:33.36 


12:17:24.24 


13:47:33.35 


12:17:24.00 


0.15 


0.24 


Mrk 463E 


13:56:02.78 


18:22:18.59 


13:56:02.80 


18:22:18.79 


-0.28 


-0.20 


IRAS 14348-1447NE 
IRAS 14348-1447SW 


14:37:38.35 
14:37:38.23 


-15:00:21.27 
-15:00:24.28 


14:37:38.38 
14:37:38.26 


-15:00:21.56 
-15:00:24.57 


-0.43 
-0.43 


0.29 
0.29 



References. — Radio Coordinates: Condon et al. 1990. CO Coordinates: This work. 



Table 4 
CO Emission Line Properties 



Source 


D L a 


2CO 


AtlFWHM 


ScoAu 




M(H 2 ) b 




(Mpc) 




(km s" 1 ) 


(Jy km s" 1 ) 


(K km s" 1 pc 2 ) 


(M S ) 


IRAS 08572+3915 


236 












IRAS 08572+3915NW 




0.0583 


250 


11.8+1.4 


1.5 x 10 9 


2.3 x 10 9 


IRAS 08572+3915SE 






(250) c 


< 4.2 d 


< 5.4 x 10 8 


< 8.2 x 10 8 


IRAS 12112+0305 


296 






39.7+4.8 


8.0 x 10 9 


1.2 x 10 10 


IRAS 12112+0305NE 




0.0728 


400 


29.8+2.4 


6.0 x 10 9 


9.0 x 10 9 


IRAS 12112+0305SW 




0.0729 


300 


9.9+2.4 


2.0 x 10 9 


3.0 x 10 9 


IRAS 13451 + 1232 


503 












IRAS 13451 + 1232NW 




0.1224 


600 


14+4 


7.8 x 10 9 


1.2 x 10 10 


IRAS 13451 + 1232SE 






(600) c 


< 12 d 


< 6.7 x 10 9 


< 1.0 x 10 10 


Mrk 463 


206 












Mrk 463E 




0.0505 


200 


6.8+0.9 


6.8 x 10 8 


1.0 x 10 9 


Mrk 463W 






(200) c 


< 2.7 d 


< 2.7 x 10 8 


< 4.0 x 10 8 


IRAS 14348-1447 


335 






40.8+4.4 


1.0 x 10 10 


1.6 x 10 10 


IRAS 14348-1447NE 




0.0823 


350 


15.6+2.2 


4.0 x 10 9 


6.0 x 10 9 


IRAS 14348-1447SW 




0.0827 


300 


25.2+2.2 


6.4 x 10 9 


9.6 x 10 9 



a Luminosity Distance, calculated assuming Hq — 75 km s 1 Mpc 1 , qo — 0.5, and A — 0. 
b Calculated assuming a. — 1.5M© [K km s^ 1 pc 2 ] -1 . 

c Adopted emission line width for upper limit determinations of Sco^v, L co , and M(H2). 
d i.e., 3Sco(rms)Ai;FWHM where 5co(rms) is the root-mean-squared flux density. 
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Table 5 
Estimation of Md yn /L[ 



Source 


-Rradio 
(PC) 


(Jy) 


/l00(im 

(Jy) 


7b b 

(K) 


Rco 

(PC) 


May" 

(M ) 


Mdyn/LcO 

(Mq [K km s" 1 pc 2 ]" 1 ) 


IRAS 08572+3915NW 


40 


7.43 


4.59 


120 


130 


1.3 x 10 9 


1.2 


IRAS 12112+0305 
IRAS 12112+0305NE 
IRAS 12112+0305SW 


65 
70 


8.50 


9.98 


65 


270 
180 


7.2 x 10 9 
2.7 x 10 9 


1.6 
1.8 


IRAS 13451+1232NW 




1.92 


2.06 


70 


240 


2.0 x 10 10 


2.5 


Mrk 463E 




2.18 


1.92 


80 


120 


1.0 x 10 9 


1.5 


IRAS 14348-1447 
IRAS 14348- 1447NE 
IRAS 14348- 1447SW 
Average 


85 
110 


6.87 


7.07 


70 


220 
300 


5.3 x 10 9 
5.5 x 10 9 


1.6 
1.0 
1.6 



a Calculatcd using Rqo- 



References. — Far infrared flux densities - Kim & Sanders (1998), except Mrk 463 (IRAS Faint Source Catalog); 
8.4 GHz radio sizes - Condon et al. (1991). 



Table 6 

Flux Ratios of Multiwavelength Emission from Component Nuclei in ULIGs 



Source 


fco 


ft. 


/Pac 


/8.4GHz 


/l.4GHz 




ratio 


ratio 


ratio 


ratio 


ratio 


IRAS 08572+3915 NW/SE 


> 2.8 


2.2 


8.7 






IRAS 12112+0305 NE/SW 


~ 3.0 


2.5 


0.7 


4.5 


4.6 


IRAS 13451+1232 NW/SE 


> 1.2 










Mrk 463 E/W 


> 2.5 


17.2 






100 


IRAS 14348-1447 SW/NE 


~ 1.6 


2.4 


1.5 


1.6 


1.7 



References. — Ha: Kim et al. (1995), except Mrk 463 (Chatzichristou & 
Vanderriest 1995) and IRAS 12112+0305 (Colina et al. 2000). Paa: Murphy 
etal. (2001). 8.4 GHz: Condon et al. (1991). 1.49 GHz: Condon et al. (1990), 
except Mrk 463 (Mazzarella et al. 1991). 
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Table 7 

Star Formation Rate Estimates and Gas Consumption Timescales 



Source 




Star Formation Rates 




Gas Consumption Timcscalcs a 




(1) 
from Lfir 
(M yr" 1 ) 


(2) 

from Z,i.4GHz 
(M yr" 1 ) 


(3)° 
from S gas 
(M yr" 1 ) 


(4) 
from Mh 2 
(M yr-f) 


(1) 

from Lfir 
(10 7 yr) 


(2) 

from Li.4GHz 
(10 7 yr) 


IRAS 08572+3915 
IRAS 08572+3915NW 
IRAS 08572+3915SE 


180 


60 


180-460 


60 

<20 


1.3 


4.1 


IRAS 12112+0305 
IRAS 12112+0305NE 
IRAS 12112+0305SW 


270 


320 
260 
60 


880-2460 
680-2110 
200-350 


300 
220 
80 


4.4 


3.4 
5.2 


IRAS 13451 + 1232 
IRAS 13451+1232NW 
IRAS 13451 + 1232SE 


180 




1080 c 


290 
<250 


6.6 




Mrk 463 
Mrk 463E 
Mrk 463W 


30 




60 c 


30 
<10 


3.1 




IRAS 14348-1447 
IRAS 14348-1447NE 
IRAS 14348-1447SW 
Average 


300 


550 
200 
350 


1130-2550 

450-980 

680-1570 


390 
150 
240 


5.2 
4.1 


2.9 
2.7 
3.7 



a Scalcs inversely with a. 

b The lower SFR is calculated using size estimate 1, and the upper SFR is calculated via size estimate 2. 
c Calculatcd using size estimate method 2 only. See §5. 



Table 8 

Estimates of Surface Densities, Column Densities, and Volume Densities' 1 



Source 


Measured 8.4 GHz size constraint 
(Size Estimate 1) 




D co oc L' co /T bb Av 
(Size Estimate 2) 




JJco n p 

(po &) (»-) m 


Deo 
(pc) 


U n 





IRAS 08572+3915NW 


80 


4.4 x 10 5 


2.9 x 10 2S 


8.6 x 10 3 


250 


4.5 x 10 4 


3.0 x 10 24 


270 


IRAS 12112+0305NE 
IRAS 12112+0305SW 


130 
180 


6.7 x 10 6 
1.2 x 10 5 


4.4 x 10 25 
7.7 x 10 24 


7.7 x 10 3 
9.7 x 10 2 


540 
360 


3.9 x 10 4 
2.9 x 10 4 


2.6 x 10 24 
1.9 x 10 24 


110 
120 


IRAS 13451+1232NW 










480 


6.6 x 10 4 


4.3 x 10 24 


210 


Mrk 463E 










230 


2.4 x 10 4 


1.6 x 10 24 


160 


IRAS 14348- 1447NE 
IRAS 14348-1447SW 


170 
220 


2.8 x 10 5 
2.7 x 10 5 


1.8 x 10 2S 
1.8 x 10 25 


2.5 x 10 3 
1.9 x 10 3 


450 
600 


3.8 x 10 4 
3.3 x 10 4 


2.5 x 10 24 
2.2 x 10 24 


130 
80 



a See Section 6.3.3 for an explanation of the computations. Note that, with exception to Dco, these quantities scale 
with a. 

References. — 8.4 GHz sizes - Condon et al. (1991). 



This figure "Evans.figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0208541vl 



This figure "Evans.fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0208541vl 



This figure "Evans.fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0208541vl 
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